Vm,n 




Vm,n 




0m, n 




4>i,i 





01,9 


-1.60 


01,17 


.64 


01,2 


-.41 


01,10 


-1.53 


01,18 


1.19 


01,3 


-.76 


-01,11 


-1.4 


01,19 


1.8 


01,4 


-1.05 


^1,12 


-1.21 


01,20 


2.47 


01,5 


-1.28 


01,13 


-.96 


01,21 


3.2 


i>i,e 


-1.45 


01,14 


-.65 


01,22 


3.99 


01,7 


-1.56 


01,15 


-.28 


01,23 


4.84 


01,8 


-1.61 


01,16 


.15 


01,24 


5.75 



Table 1: Dimensions of primary operators in the (3,25) model. 



(p,q) 





$ 


(p,q) 





$ 


(2,21) 


01,4 


01,1 ■ 


■ -01,10 


(7,62) 


01,13 


01,1 ■ 


■ ■ 01,20 


(3,25) 


01,11 


01,1 ■ 


■ -01,19 






02,7- 


■ ■ 02,28 


(3,26) 


01,5 


01,1 ■ 


■ ■ 01,20 






03,16 


■ ■ ■ 03,37 


(5,32) 


01,9 


01,1 ■ 


■ -01,14 


(8,47) 


03,17 


01,1- 


■ -01,14 






02,5- 


■ ■ 02,21 






02,4- 


■ ■ 02,20 


(6,37) 


01,8 


01,1 ■ 


■ ■ 01,14 






03,10 


■ ■ ■ 03,26 






02,4- 


■ ■ 02,20 






04,16 


■ ■ ■ 04,23 






03,11 


■ ■ ■ 03,18 


(8,67) 


03,28 


01,1 ■ 


■ ■ 01,19 


(6,55) 


01,14 


01,1 ■ 


■ ■ 01,20 






02,6 ■ 


■ ■ 02,27 






02,8- 


■ ■ 02,29 






03,14 


■ ■ ■ 03,36 






03,17 


■ ■ ■ 03,27 






04,23 


■ ■ ■ 04,33 



Table 2: Allowed perturbations of the first nine minimal models. 
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So far, our treatment has been completely general. All of the above analysis applies 
to any minimal model solution of conformal turbulence. To proceed further we need 
to examine specific models. 

Let us briefly summarise what we want to find. Starting from a minimal model 
(p, q) satisfying Polyakov's constraints (5) and (7) or (8), with a certain 'stream func- 
tion' operator ?/>, we choose different perturbing primary operators $ from the model. 
Denoting by ^ the minimal dimension operator appearing in the OPE we look 

for the $'s that will satisfy condition (25). To do this we rearrange (25) to give 
A$ < 1 + — 2A^. Let us now define k = 1 + A^ — 2A^. We then calculate the 
value of k for different choices of We have done this for the first few models for 
both the constant enstrophy and the constant energy case. In general, the value of k 
depends not only on the model but also on the perturbation operator. Although it is 
not constant for a given model it varies within a well-defined narrow range. We can 
then easily identify those operators with dimension less than k. For example the value 
of k for the (3,25) model is either 2.19 or 2.2. We find that the fields ?pi,ii ■ ■ ■ , "01,19 
have dimension less than this (see Table 1). We have investigated the first few low- 
lying models. The detailed results are shown in Table 2. The value of k lies roughly 
between 1 and 3. Hence we can conclude that, at least for the models investigated, all 
the solutions of (28) also satisfy (25), an example being the case $ = ip mentioned 
above. The converse statement is not necessarily true. 

In conclusion, it would certainly be worthwhile extending the approach in this 
paper to find all possible perturbations that leave the inviscid Hopf equation invariant, 
at least to leading order in perturbation theory. Since such perturbations also break 
conformal symmetry, one can view our results as widening the possible solutions of two- 
dimensional turbulence. When conformal symmetry is broken in this way, one natural 
question is whether the parameters in the models (e.g. A , the scaling dimensions, or 
even the viscosity v ) undergo renormalization group flow, and if so what might be the 
fixed points? A related outstanding problem is to give a physical interpretation to the 
various perturbing operators discussed in this paper. The results presented in [6] could 
provide a starting point for such an investigation. Such an investigation could place 
additional constraints on the allowed set of perturbing operators listed in Table 2, for 
example that A$ > 0. 
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for all solutions except (2,21), for which the lower limit is ^y. For the constant energy 
case, the corresponding range is 

— < R < 3-V8 (33) 

22 v ; 

We have computed the values of Ao for these limits, and the corresponding ranges are 

given by, respectively, 

- 2.025 < A < -1.5 (34) 
(with the lower limit as -2.14 for the (2,21) model) and 

-1.37 < A < -1 (35) 

Now we go back to equation (28). For negative A$ the second term is positive. If for 
a moment we take tp to be the absolute minimal operator in the model then we note 
that (28) will automatically be satisfied if we can fulfill the stronger constraint 

A + 2 > (36) 

We note that this is manifestly satisfied by all the solutions of the constant energy 
constraint. For the constant enstrophy case, the only possible violations arise from 
the (2,21) model, where Aq = 2.14, and those models with —2.025 < Aq < —2, ie 
0.1 < R < 0.101. In each of these, if ip is any operator with dimension < —2, then 
(28) will be violated if there is any perturbing operator $ which contains ip in its OPE 
$ $ such that 

A$ > 1 + ^A v (37) 

Explicit calculation for the (2,21) model shows that this is not the case. In fact, this 
last condition seems to be very unlikely to be met by any minimal model. All the 
models which we have tested seem to indicate that (28) is satisfied anyway, regardless 
of whether Aq < —2. Unfortunately we have not as yet been able to come up with a 
rigorous proof to that effect. 

To summarise, we have found that all primary operators with negative dimension 
and some of those with positive dimension < 1 can be used to perturb a general non- 
unitary CFT without introducing ultraviolet divergences. 

With regard to the solutions of conformal turbulence, we now return to equation 
(25) and try and verify if there is at least a subset of these operators that also leaves 
the Hopf equations invariant to 0(A 2 ). 
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an overall negative result, violating (28). Note that this latter result is quite general 
and in fact applies to any CFT, not just to those satisfying conformal turbulence. The 
remaining region < A$ < 1 is more delicate to analyse. It appears that the result is 
dependent on the particular model under consideration and that no general statement 
can be made here. We illustrate with the example of the (3,25) model. In this model, 
we have two such operators: ^1,16 an d ^1,17, with dimensions 0.15 and 0.64 respectively 
(see Table 1). First let $ = ^1,16- The minimal dimension field in its OPE is ?/>i g, 
with dimension A v = —1.6. Evaluating the lhs of (28) we find that the inequality is 
indeed satisfied. Next we choose $ = Vi,i7- This again gives tp = Vi,9- Inserting 
the values for the corresponding dimensions, we see that (28) is now violated. Thus 
we obtain opposite results for these two operators. 

We next consider the final case where A$ < 0. Consider the minimal dimension Aq 
appearing in the non-unitary minimal (p, q) model. This is given by 

Defining R = p/q and rearranging equation (29) we obtain 

A = i[2- J R-i(l-l/ 9 2 )] (30) 

Now the smallest value of q occuring in the solutions of conformal turbulence is 21, 
making the 1/q 2 term negligible. Thus, 

A ~ \ [2 - R - 1/ R] (31) 

The percentage error introduced by using the latter expression is less than 0.3% for 
the (2,21) model and decreases rapidly with increasing q for subsequent models. Since 
p < q, the physical domain of R is (0,1). In this range Ao is a monotonically increasing 
function as illustrated in fig 1. For large q, the limit R = 1 corresponds to unitary 
minimal models, for which q = p + 1. There we see that Ao = 0, ie the minimal 
dimension operator in such models is the unit operator, with conformal dimension 0. 
In other words, no negative dimensions appear in such theories. As far as the minimal 
models of conformal turbulence are concerned, we saw in the introduction that they 
are always of necessity non-unitary. For these, Lowe [2] derived some inequalities that 
R should satisfy. For the constant enstrophy case, these are 

< R < 4 -a/15 (32) 
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As emphasised in [7], one has to consider the possibility that new ultraviolet di- 
vergences can arise in correlation functions of the perturbing operator. To see how to 
avoid such divergences, consider correlation functions of the perturbing operator $ in 
the perturbed theory. Similar to (13) we obtain 

< ^(z 1 ,z 1 )^(z 2 ,z 2 )-- - >s = < ^(z 1 ,z 1 )^(z 2 ,z 2 )-- - > s . 
+A J <$(z,z)$(z 1 ,z 1 )$(z 2 ,z 2 )---> s . d 2 z +0(\ 2 ) (26) 

The ultraviolet behaviour of the integrand in the first order term is determined by the 
operator product expansion of $ with itself 

$ $ ~ | fl |2(A v -2A«) ^ ^7) 

where ip is the minimal dimension operator in the OPE in equation (27). There will 
be no new ultraviolet divergences provided 

A v - 2A,e + 2 > (28) 

Now, in a unitary theory, since all the conformal dimensions are positive, it is clear 
that this condition will automatically be satisfied if A$ < 1, as pointed out in [7] i.e. $ 
is a so called relevant perturbation. However, in non-unitary theories the situation is 
further complicated by the presence of fields with negative dimensions, and more care 
is required. Our immediate task is to find general solutions of (28). Later we also look 
for sub-classes of these that satisfy (25), thus leaving the Hopf equation invariant to 
0(A 2 ). 

For illustration, we first consider some special cases. For example, if we perturb 
with the stream function operator ip itself we have $ = ip and ip = <p. Then 
we can make use of either (7) or (8) to eliminate A v from (28). The result is A^ < 
— | for the constant enstrophy case and A^ < for the constant energy case. We 
recall that the corresponding conditions imposed on ip by conformal turbulence are 
A^ < —1, and A^ < — | in each of these cases respectively. Hence, we can perturb 
any solution of conformal turbulence with the stream function operator ip without 
introducing ultraviolet divergences in the 2-point functions of ip. 

Let us now be more general. We consider separate cases. For positive A$, we note 
that, since the unit operator always occurs in the OPE $ the maximum possible 
value of the first term A v is 0. Thus, if A$ > 1 the second term will dominate and yield 
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the action of the final d z are of the form bDip^ i. Performing the integrations, we finally 
obtain for the integral / 

/ = (2tt) 2 lim A T l&^A^-A^-Ao + l) l^tA.rA,,-^) p ^(^j) ( 2 1) 
a— >-0,6 — >-0 . . . , 

To avoid ambiguities in the ordering of the limits a — > and ft — >■ 0, we shall take both 
limits simultaneously. The effect of taking the simultaneous limits is that a and b can 
be considered to be essentially the same cut-off parameter, and therefore we can write 

/ ~ lim A^|a| 2 ( A *.'- A *- 2A * +1 ) D VvM) (22) 

k,l 

Note that the dimensions of the intermediate fields ipij drop from the final result. The 
exponent of \a\ is crucial. To avoid introducing singularities the perturbing operator $ 
has to be such that none of the fields ipk,i gives rise to a negative exponent. In other 
words, for any sensible perturbation we have the necessary condition 

A*,/ - A* - 2A^ + 1 > (23) 

The limit a — ¥ means that we need only keep the term with the smallest exponent 
of \a\. Let us denote the corresponding lowest dimension field ijjhj by ^ . We can then 
write 

/ rsu | a |2(A*-A*-2A^ + l) D y (24) 

This is the direct analogue of equation (4). Here /, as defined by (14), is u dressed by 
the conformal perturbation . D is a (2,2) parity even operator made up from L-i, 
L_2 and Z_2, and ^ is the minimal dimension operator appearing in the OPE iptp^i 
analogous to <f> of equation (4). 

We thus find that, in the limit a — Y 0, / will vanish for perturbations such that 

A* - A* - 2A^ + 1 > (25) 

This is analogous to condition (5). In that case, we conclude from equation (13) that 
the Hopf equation will still be satisfied, at least to first order in A. Hence, the minimal 
model solutions of conformal turbulence are stable to such perturbations. In other 
words, although perturbing these CFT's generally breaks their conformal invariance (ie 
they are no longer CFT's) they still remain solutions of the Navier-Stokes equations. 
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where — in lining denotes spatial averaging over the cutoff a, i.e. if we express a = 
\a\e tda then one must integrate over 6 a before taking the limit a — Y Substituting this, 
we then obtain 

/ = lim / d 2 uj 0(z, z, z\ z ) ip{z, z)ip(z\ z')<$>(lj, uj) , (17) 

where 

0(z,z,z',z') = d z ,d 2 d z - d- z ,d 2 d- z 

and z' = z + a, z' = z + a. As it stands this integral is ultraviolet divergent and needs to 
be regularised. Indeed it is because of this that we expect to obtain generally different 
kinds of operators than the one found by Polyakov in equation (4) . Following [7], we 
do this by inserting step functions Q[(z — oo)(z — uj) — b 2 ] and Q[(z' — uj)(z' — Co) — b 2 ] 
into the integral, (and again taking spatial averaging over the cutoff b ). Next we pull 
out a d z and act with it on the first G-function. This gives a ^-function, which enables 
us to perform the radial integration. The result is 

h = lim r de a r de b \b\ d z (d z ,d- z - d z ,d z ) v^^MW^)] , (is) 

a-S-0 JO JO 

where b = \b\e tdb , If, is the regularized integral / and the square brackets mean that we 
have to take the OPE of the last two operators first. Performing this OPE, we obtain 

il>(z,z)*(u,u) - J2J2\b\ 2{A ^- A ^Az,m (19) 

i j 

The limit a — Y then forces us to take the next OPE. Doing this we find 

i>{z'^)[i>{z,z)^{u,Q)] ~ Y. |6| 2 ( A -- A ^)|a| 2 ( A ^- A -- A ^ [Mz,z)] (20) 

where, following usual convention, [ipk,i] denotes the conformal class of ipk,i- The de- 
scendants of the primary fields are crucial. As in Polyakov's treatment, we need to 
go beyond the leading term to obtain a non-vanishing result. It is clear from parity 
considerations that the effect of the operator (d z jd z — d z >d z ) in (18) is to yield a 
lowest level term of the form (ab — ba) Dip/,,1, where D is a dimension (2,2) parity even 
operator such as L_2-£-2, L^L 2 ^, L_2L 2 _ 1 etc. Equation (19) implies that, due to the 
angular integration, a non-vanishing result will be obtained only if we can cancel the 
factor of e t6b . This means that we need an extra b. Hence the only terms that survive 
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2 Perturbing the solutions of conformal turbulence 



We now wish to investigate the effect of perturbing such a (non-unitary) CFT by some 
primary operator in the theory. In particular we are interested in finding out whether 
the perturbed theory still satisfies the Hopf equations. This will give some indication 
on the 'stability' of these solutions. 

We proceed in a similar manner described by Cardy [7]. Consider the effect on the 
'action' describing the CFT of perturbing by a certain primary operator $ with scaling 
dimensions (h, h): 

S = S*-\ J $(z,z)d 2 z , (12) 

where S* is the fixed point action and A is a coupling constant with conformal di- 
mensions (1 — h, 1 — h). Evaluating the correlation functions in the Hopf equation (2) 
perturbatively in A, we obtain (in complex notation) 



< u(z 1 ,z 1 )u(z 2 , z 2 ) ■ ■ ■ > s = < u(z 1 ,z 1 )u(z 2 ,z 2 ) ■■■>$• 
+A J <<f>(z,z)u(z 1 ,z 1 )u(z2,Z2)--- >s* d 2 z + 0(\ 2 ) (13) 

The first term is simply the unperturbed Polyakov correlator evaluated in the minimal 
model. Thus, by (4), it vanishes in the usual Polyakov ansatz, namely equation (5). 
Assuming this condition, we want to determine the effect of the second term. It is clear 
that the only contributions which do not obviously vanish as a result of equation (5), 
are those obtained by evaluating the integral 

I = \J ®{z,z) u d 2 z (14) 

We can use equation (3) to substitute for u, obtaining 

I = -A J d 2 uj [e a/3 d a 'ip(z,z)df ) d 2 'iij(z,z)} , (15) 

Equation (15) represents the dressing of the inertial term by the perturbing operator, 
and as before, the term in brackets is to be understood in terms of its point-splitted 
definition 

e a pd a ip(z, z)d f) d 2 tp(z, z) = lime Q/3 d Q V(^ + a , z + a)dpd 2 i)(z, z) (16) 
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One of the most important physical predictions of conformal turbulence is the 
exponent of the energy spectrum. This is found to be given by 



E(k) ~ k 4A ^ +1 (9) 

Now, since the dimension of ijj is model-dependent, the predicted exponent is not 
unique. The experimental value lies between 3 and 4, but there are many models 
which give this result. Some additional constraints are needed to identify the correct 
theory of conformal turbulence. This has not been done so far. 

The simplest solutions are minimal models. In general they will always be non- 
unitary. This means that there will be some primary operators with negative dimen- 
sions in the theory. To see this we note that at least the stream function operator 
ip will always have negative dimension. This is because equations (5) and (7) taken 
together imply that A^ < — 1 while (5) with (8) give A^ < — |. In unitary models the 
unit operator, with dimension zero, is always the lowest dimension operator. However, 
for the non-unitary models, due to the negative dimensions, this is no longer the case. 
Some pertinent facts about minimal models are summarised here. 

The (p, q) minimal model (where p and q are co-prime positive integers, with p < q) 
contains |(p — l)(q — 1) degenerate primary operators ip m ,n, where 1 < m < p and 
1 < n < q, with conformal dimensions given by 

A mn = {pn - qm)2 - (p - q)2 (10) 
Apq V ' 

These primary operators satisfy the following fusion rules: 



i j 

where % runs from \m\ — m2\ + 1 to Min(mi + — l,2p — m\ — — 1), and is 
odd if mi + ™<2 is even and vice- versa, and similarly for j. The D\ h ^ ■. are the 

structure constants of the CFT. The symbol [ipij] denotes the conformal family of ipij. 
The operators ip ryS and ipp-r, q - s have the same dimensions, and they are identified. It 
is found that an infinite number of these minimal models satisfy the constraints (5) 
and (7) or (8), the simplest being the (2,21) model with ip = V'1,4- Some results are 
listed in the paper by Lowe [2]. 
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range: a < I. Here / is the infrared cutoff and a is the ultraviolet cutoff, which 

is determined by the viscosity. Specifically, the stream function ip is assumed to cor- 
respond to a certain primary operator of some CFT. In the limit v — Y (ie a — Y 0), 
equation (1) reduces to 

^ = -e a pd a 4>dpd' 2 'il) (3) 

The right hand side is not well-defined and needs to be regularised by point-splitting. 
The result is 

6j ~ \a\ 2A *- 4A 4(L-2L 2 _ ! - L\L_ 2 )4> (4) 

where <f> is the minimal dimension operator which appears in the OPE of ijj with itself. 
The inviscid Hopf equation will only be satisfied if to vanishes. From the above result we 
see that there are two possibilities. The first is that the operator (L_2-£ 2 _ 1 — L 2 _ 1 L-2)4> 
is identically zero. This condition is trivially fulfilled by requiring the corresponding 
CFT to be degenerate on level two. The simplest such solution is the (2,5) minimal 
model. The second possibility is that 

> 2A^ (5) 

Thus, any theory with a 'positive defect of dimensions' will solve the Hopf equations. 
A further constraint introduced by Polyakov [1] is the constant enstrophy flux condi- 
tion, 

< 6j(r)oj(0) >= const (6) 

which leads to the condition 

+ = -3 (7) 

Physically, constant enstrophy flux guarantees that conformal field theory correctly 
describes the enstrophy input at large scales being dissipated at small scales. Polyakov 
found that these requirements are fulfilled by the (2,21) minimal model. Subsequently, 
Lowe [2] and others [3-5] have shown that there are an infinite number of minimal 
model solutions of these constraints, the (2,21) model being the simplest one. Lowe 
also considers an alternative constraint, the constant energy flux condition, which yields 

A^ + A^ = -2 (8) 

This also has an infinite number of solutions. The constant energy flux condition 
allows for the possibility of an inertial range in which energy cascades from small to 
large scales as envisaged by Kraichnan [10], Leith [11] and Batchelor [12]. 
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1 Introduction 



Following the recent proposal by Polyakov [1], that turbulent flow in two spatial di- 
mensions may be understood in terms of certain non-unitary conformal field theories 
(CFT) (at least in the inviscid limit), there has been renewed interest in such systems. 
It was quickly realized that Polyakov's results do not single out a unique CFT, rather 
there appears to be an infinite number of suitable theories [2-5], which even includes 
the case when boundaries are present that lead to non-vanishing one-point functions 
[1]. Although the appearance of an infinite number of solutions might be an unde- 
sirable feature, the power and elegance of CFTs in this context have led to further 
investigations and predictions concerning two dimensional turbulence [6]. 

In this paper we shall consider the situation when any of the non-unitary CFT 
solutions to 2-dimensional turbulence is perturbed by some primary field of that CFT. 
Whilst the consequences (see e.g. [7]) of perturbing unitary CFT's in this way have been 
well documented (perhaps the most notable example being the so called C-theorem of 
Zamolodchikov [8]), less is known concerning non-unitary CFT's [9]. Furthermore we 
shall investigate the effects that such perturbations have on the Hopf equations that 
describe the statistical properties of the inviscid turbulent flow. In particular we show, 
by considering the first few CFT solutions to turbulence, that there always exists 
primary fields that leave the Hopf equations invariant. 

We begin by reminding the reader of some of the basic formulae underlying Polyakov's 
solution to 2-dimensional turbulence. The turbulent flow of a fluid with viscosity v is 
governed by the Navier-Stokes equations, which in two spatial dimensions take the 
form 



where tp is the stream function, related to the velocity v by v a = e a pdfii> an d to the 
vorticity u by uj = d 2 ij) . In the statistical formulation of turbulence one consid- 
ers correlation functions of relevant quantities (velocities and vorticities) . Demand- 
ing that we have a stationary probability distribution implies that the correlators 
< oj(xi) ■ ■ ■ oj(x n ) > are time-independent. This leads to the inviscid Hopf equation: 



Polyakov [1] conjectured that we could interpret these correlators in terms of an ef- 
fective conformal field theory if all the coordinates Xi are within the so-called inertial 




< u(x 1 )u(x2) •••> + < u(x 1 )u(x2) •••>-! = 



(2) 
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